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Abstract 

We construct a realistic Supersymmetric Grand Unified Theory of Flavour based on 
P SX2(7) x 5*0(10), where the quarks and leptons in the 16 of 50(10) are assigned 
to the complex triplet representation of PSL^iJ), while the flavons are assigned 
to a combination of sextets and anti-triplets of PSL2(7). Using a D-term vacuum 
alignment mechanism, we require the flavon sextets of PSL2{7) to be aligned along 
the 3-3 direction leading to the third family Yukawa couplings, while the flavon 
anti-triplets describe the remaining Yukawa couplings. Other sextets are aligned 
along the neutrino flavour symmetry preserving directions leading to tri-bimaximal 
neutrino mixing via a type II see-saw mechanism, with predictions for neutrinoless 
double beta decay and cosmology. 
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1 Introduction 



The discovery of neutrino mass and approximately tri-bimaximal (TB) lepton mixing pQ 
suggests some kind of a non-Abelian discrete family symmetry Gf might be at work, at 
least in the lepton sector. In the neutrino flavour basis (i.e. diagonal charged lepton 
mass basis), it has been shown that the TB neutrino mass matrix is invariant under 5, U 
transformations, M^ B = SM£ B S T = UM^ B U T [2|. A very straightforward argument 
shows that this neutrino flavour symmetry group has only four elements corresponding to 
Klein's four-group Zif x Z^- By contrast the diagonal charged lepton mass matrix (in this 
basis) satisfies a diagonal phase symmetry T. The matrices S, T, U form the generators 
of the group 5 4 in the triplet representation, while the A A subgroup is generated by 5, T. 
The observed neutrino flavour symmetry corresponding to the two generators 5, U may 
arise either directly or indirectly from a range of discrete symmetry groups [3|. Examples 
of the direct approach, in which one or more generators of the discrete family symmetry 
appears in the neutrino flavour group, are typically based on S4 [I] or a related group such 
as A 4 [5l[6] or PSL 2 (7) [2]. Models of the indirect kind, in which the neutrino flavour 
symmetry arises accidentally, include also A A [7J and S4 [S] as well as A 2 7 [S] and the 
continuous flavour symmetries like, e.g., 5*0(3) JTU] or SU(3) [UJ which accommodate 
the discrete groups above as subgroups [12]. For an incomplete list of models with family 
symmetries see [T3Tf2"T] and the reviews [28|[29j. 

A desirable feature of a complete model of quark and lepton masses and mixing angles 
is that it should be consistent with an underlying Grand Unified Theory (GUT) structure, 
either at the field theory level or at the level of the superstring. The most ambitious models 
which have been built to achieve this are based on an underlying SO (10) structure. This 
is very constraining because it requires that all the 16 spinor components of a single 
family should have the same family charge, comprising the left-handed fermions tp and 
the CP conjugates of the right-handed fermions ip c , including the right-handed neutrino. 
Although it seems somewhat ad hoc to add right-handed neutrino singlets to 577(5), once 
they are present it is straightforward to construct models of lepton masses and TB mixing 
that are consistent with quark masses and mixing, and there are several successful models 
of this kind for example based on A4 x 577(5) |23j . 

Despite the theoretical attractiveness of 50(10), there are very few 50(10) models 
capable of enforcing TB mixing by means of a family symmetry. It is desirable that such 
models contain complex triplet representations of the family symmetry and examples of 
such models based on the Pati-Salam subgroup of 50(10) have been constructed where 
the family group is SU(3)f [11] or A27 [9], with Yukawa couplings arising from operators 
of the form fti/j^ipj. If the triplet representations were taken to be real rather than 
complex then this would allow an undesirable alternative contraction of the indices in 
the Yukawa operator namely fyipffi 'ft H leading to a Yukawa matrix proportional to the 
unit matrix which would tend to destroy any hierarchies in the Yukawa matrix. For real 
representations, there is no symmetry at the effective operator level that could forbid such 
a trivial contraction of the two triplet fermion fields. However in principle it is possible 
to appeal to the details of the underlying theory in order to forbid the trivial contraction. 
This involves a discussion of the heavy messenger states whose exchange generates the 
operator, for example as was done recently in the 54 x 50(10) model in [8]. 
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A general problem with all 5*0(10) models in which the quarks and leptons form 
triplet representations of the family symmetry (necessary to account for TB mixing) is 
that the top quark Yukawa coupling arises from a double flavon suppressed operator of the 
form ^i/j^ipj. The situation improves somewhat if one considers flavons also in sextet 
representations, as pointed out in [2]. For example, introducing a two index anti-sextet 
of SU(3)f, x %1 ') the lowest order Yukawa operators become, x^ip^. If the anti-sextet 
flavon x has a VEV aligned along the 3-3 direction = VS^Sj^, then this operator 
with a coefficient y/M would imply a third family Yukawa coupling of yV/M with the top 
quark Yukawa coupling of 0.5 implying V/M « 0.5y~ l , which has acceptable convergence 
properties. 

Flavon sextets are therefore well motivated as the origin of the third family Yukawa 
couplings in SO (10) models. Although this is possible in the case where the family group 
is SU(3)f [UJ, it is not possible for the discrete groups A 27 [9] or S4 jS] for the simple 
reason that these groups do not admit sextet representations. The smallest simple discrete 
group which contains complex triplets and sextet representations is PSLi2(7), which is the 
projective special linear group of two dimensional matrices over the finite Galois field of 
seven elements. PSL 2 (7) contains 168 elements and is sometimes written as £(168) [3D] . 
The relationship of PSL 2 (7) to some other family symmetries that have been used in the 
literature is discussed in [26" | l3~T j l3"2] . 

In a recent paper [2] we developed the representation theory of PSL 2 (7) for triplets and 
sextets in a convenient basis suitable for applications of PSL 2 (7) as a family symmetry 
capable of describing quark and lepton masses and mixing angles in the framework of 
5*0(10) models. We showed how the triplet representation given in terms of the standard 
generators A, B in [21] may be related to four PSL 2 (7) generators S,T,U,V. In such 
a basis the subgroup structure PSL 2 (7) D 5*4 D A 4 just corresponds to the respective 
generators being 5, T,U,V D S,T,U D 5, T. 

The purpose of the present paper is to construct a realistic Supersymmetric Grand 
Unified Theory of Flavour based on PSL 2 (7) x 5O(10). We require the sextets to be 
aligned along the 3-3 direction to account for the large third family Yukawa couplings, 
while we shall make use of anti-triplet flavons, whose VEVs are aligned along the columns 
of the TB mixing matrix, to account for the first and second family quark and lepton 
masses and mixings. It turns out that in PSL 2 (7) it is easier to obtain the sextet flavon 
vacuum alignments first using the D-term approach to vacuum alignment discussed in [3] , 
but here applied to flavon sextets rather than flavon triplets. In this way we can obtain 
sextet vacuum alignments along the 3-3 direction, suitable for the top quark Yukawa 
coupling, if we assume two relations (possibly arising by virtue of a higher symmetry) 
amongst different quartic sextet combinations appearing in the flavon potential. Sextet 
vacuum alignments along the 5, U preserving directions, suitable for reproducing the 
neutrino flavour symmetry of TB mixing in a direct way as a subgroup of PSL 2 (7), 
can be naturally obtained. Once the sextet flavons have been aligned, the anti-triplet 
flavons are then aligned against the pre-aligned sextet flavons, in particular using the 
5, U preserving sextet flavons, leading to triplet flavon alignments along the columns of 
the TB mixing matrix as mentioned. 

In general the neutrino mass and mixing in this model can arise from either the 
type I see-saw or the type II see-saw or both. The discussion of the type I see-saw 
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Figure 1: Scatter plots for m ee and \rrii\ against m m \ n for the type II see-saw model based 
on PSL 2 (7) x 50(10). 

approach follows along the lines of the models in [9l[ll] based on constrained sequential 
dominance [33], leading to the indirect type of models, since the triplet flavon aligned 
along the third column of the TB mixing matrix breaks both S and U . Here we shall 
focus on the new possibility offered by the PSL 2 (7) x S'O(IO) model of the type II see- 
saw mechanism where the sextet flavons aligned along the 5*, U preserving directions enter 
the neutrino sector, thereby preserving these generators leading to the neutrino flavour 
symmetry being reproduced in a direct way. As a preview of our results, we shall find that 
the type II see-saw model leads to TB mixing in the neutrino sector with a mass spectrum 
spanning the range from hierarchical or inverse hierarchical, up to the quasi-degenerate 
region. The resulting predictions for neutrinoless double beta decay mass parameter m ee 
and the total sum of physical neutrino masses JY \ m i\ relevant for cosmological hot dark 
matter are both shown in Fig. [T] as a function of the lightest physical neutrino mass m m i n 
using a double logarithmic scale. These results are in sharp contrast to the case of the 
type I see-saw possibility with constrained sequential dominance which would lead to a 
strong neutrino mass hierarchy with m m [ n = \mi\ ~ 0. Later on we shall analyse the type 
II see-saw results in more detail using linear scales. 

The layout of the remainder of the paper is as follows. In Section 2 we present the 
symmetries and superfield content of the SUSY PSL 2 (7) x S'O(IO) model and discuss 
the desired vacuum alignments and the resulting fermion mass matrices. In Section 3 
we discuss the type II neutrino phenomenology in the model. Section 4 is devoted to a 
detailed analysis of the D-term vacuum alignment for the sextet and ant-triplet flavons. 
Section 5 concludes the paper. 
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2 Fermion masses in the PSL 2 (7) X £0(10) model 



2.1 The desired vacuum alignments 

As discussed in the Introduction, the sextet flavons are used in two ways in the model. 
The sextet xlip aligned along the 3-3 direction is responsible for the third family Yukawa 
couplings, including that of the top quark, via operators like (dropping the Higgs fields) 

vwiUi • (2.i) 

The sextets Xtb aligned along the S and U preserving directions are responsible for the 
effective Majorana couplings, via operators like (again dropping the Higgs fields) 



J- 



(2.2) 



where the first term above contributes to the type II see-saw mass while the second term 
is responsible for the heavy right-handed neutrino masses. As discussed in [2], the sextet 
flavon 3 by 3 matrices x which enter the above couplings to quarks and leptons, are related 
to the six component column vector 



X = (Xi,X2,X3,X4,X5,XG) rj 



(2.3) 



by, 
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where we adopt the notation of the 'Atlas of finite groups" 

b 7 = h-l-iV7) . 



(2.5) 



Mass matrices which are of the tri-bimaximal form are obtained from alignments of the 
form 

(Xtb) = (0 , , , a 4 , a 5 , a e ) T . (2.6) 

The fully realistic type II model will require three flavons of this type which we label as 
Xtbi w here p = 0,1, 2. Explicitly their alignments read 
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Table 1: The particle content of an SO (10) model with the family symmetry PSL2(7). 

In order to generate a Yukawa matrix which gives mass to only the third generation (top 
quark, bottom quark or tau lepton), we need the alignment 

(Xtop) « • (1 , *s/3 , - iV3/b 7 , - V3/b 7 , - V2 , Of . (2.10) 

To obtain the alignments of Eqs. ( l2.7H2.10j) it is necessary to study the PSL 2 {7) symmetric 
potential for the flavon sextet, which we postpone to Section HI 

In addition the model also relies on the anti-triplet flavon fields 023 and 0123 whose 
VEVs become aligned along the directions 

, 23 > oc -L . | 1 J , (0 123 ) oc ■ I 1 I , (2.11) 

where again we postpone discussion of these alignments until later. 

2.2 The symmetries and operators of the model 

Having indicated the desired vacuum alignment of the flavon fields Xto P > Xtb? 023; 0i23> 
we can now formulate an 5*0(10) model of fermion masses using the P5L 2 (7) family 
symmetry. Table [1] lists the particle content together with all transformation properties. 
Here tp denotes the 16 of 5O(10) which incorporates the SM fermions. Even though 
our model is based on 5O(10), it is convenient to distinguish the left-handed and the 
right-handed components by showing the decomposition into Pati-Salam representations. 
Using the order SU(4) x SU(2) L x SU(2) R we can write 

16 -> (4,2, 1) + QA^l • (2-12) 

H w and Hj%? are the 5O(10) Higgs fields whose SU{2) L doublet components enter the 
Yukawa couplings. In Pati-Salam language the relevant components that acquire an 
electroweak VEV are 

H 10 (1,2,2) , H m -> (15,2,2) . (2.13) 
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Figure 2: The diagrams which generate the Yukawa couplings. In Pati-Salam language, 5J7(4) x 
SU(2)l x SU(2)r, the relevant components of the SO(10) representations are: ip — > (4,2,1), 
r -> (4, 1, 2), £T 10 -> (1, 2, 2),H m ^ (15, 2, 2). 



The latter gives rise to the Georgi-Jarlskog (GJ) factor of 3 in the (2,2) entry of the lepton 
mass matrix [291135] . With the above specified fields, the leading Yukawa operators take 
the form, 
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(2.14) 



where we have only given the PSL2(7) indices The corresponding diagrams are 

depicted in Fig. [2j Note that we have inserted a PSL2(7) singlet flavon £ in order to 
additionally suppress the third term. Inserting the flavon VEVs, the first term (a) fills in 
the (3,3) entry of the Yukawa matrix, the second (b) generates non-vanishing entries in 
the 2-3 block, and the third (c) enters everywhere in the Yukawa matrix except for the 
(1,1) component. Including sextet and anti-triplet flavons, we find the following structure 
of the Yukawa matrices, 
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H w . (2.15) 



To identify the Yukawa matrices for the different sectors Y u,d,e ' u we note that the Higgs rep- 
resentations if io and Hj^ each have two Higgs doublets in them, according to Eq. (12.131) . 
The two MSSM doublets H u , originate below the GUT scale, from one linear combi- 
nation of these up-type doublets and one linear combination of the down-type ones which 
remain almost massless. The orthogonal linear combinations acquire GUT scale masses 
just as the colour triplets and other non-MSSM states. Electroweak symmetry is broken 
after the light MSSM doublets H u , H d acquire VEVs v u ^ and they then generate the 
fermion masses. Since the dominant contribution to the 2-3 block of the Yukawa matri- 
ces arises from the components of H u , Hd coming from Hj^, these entries will receive a 
relative Clebsch factor of 3 for the leptons as compared to the quarks. Thus, ignoring 
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messenger effects for the moment, the resulting Dirac mass matrices take the form 

~ [ I with ( fl = 1 _ *> ru ' d ' (2.16) 

1 e s _ ae 2 l J M H { a = -3 for e, v , 

where, for simplicity, we have chosen 

KXtop)! n , I (023) I 2 2 -1(0123) 1(0 3 ,„ , 

Notice the zero in the (1,1) entry whose presence allows to accommodate the phenomeno- 
logically successful Gatto-Sartori-Tonin relation |36j. Moreover, the up and the down 
quark sectors can have independent messenger masses M -> M u , M d , where M u « 3M d 
so that two different expansion parameters e — > e u ,€d with e u ss ed/3 are introduced as 
in [TT] . Numerically we need e M ~ 0.05 and ~ 0.15. See [37] for a detailed discussion of 
the numerics including \ 2 fits. Thus the following mass matrix structures for the quarks 
are quite elegantly achieved: 

™ u ^ I c 3 ,2 j | KXtop)| d I 3 2 ,2 | l(Xtop)| roifi^ 

m l u el ~l u J ~mT Vu ' m [y 4 \ 7 rf J ~mr Vd ' (2 - 18) 

with similar considerations in the lepton sector leading to m e,y being the same as m d,u 
apart from the GJ factors of 3 in the 2-3 block, as in Eq. f 1 2 . X 6 j) . The (left-handed) quark 
mixing angles in the up and the down sector are then calculated as Ofy ~ e U) 4, ~ e^ d , 
and 6±f ~ d , so that the down sector gives the dominant contribution to a viable CKM 
matrix. Regarding the leptons, the mixing angles are changed by G J factors of 3 as follows 
@i2 ~ e <i,u/3, 023 ~ 3e| n , and 9^ ~ e\ u . The effect of the charged lepton mixing angles 
on the PMNS matrix is small so that tri-bimaximal mixing is only perturbed within the 
experimentally allowed range [38J. 

In this model the successful mass matrices in Eq. (12.181) are achieved in a very natural 
way since, with the inclusion of the singlet flavon £, the first row and column is cubic 
in the messenger mass, while the 2-3 block is quadratic and the 3-3 element involves the 
Higgs messenger mass Mu and so is universal. 

Turning to the Majorana sector different components of the field are responsible 
for the masses of the left-handed and the right-handed neutrinos. Analogous to the 
convention used to distinguish the left-handed from the right-handed doublet components 
of if), we introduce a notation that allows us to tell apart the left-handed and the right- 
handed triplet components of Aj^, i.e. 



126 -> (10,3,1) + (10,1,3) + ••• , (2.19) 




where the ellipsis denotes the rest of the decomposition of the 126 of 5O(10) into rep- 
resentations of 5C/(4) x SU(2)l x SU(2) r . With these remarks, the leading Majorana 
operators take the form 

£ Maj ~ i fl tydf^ATx + ^x^rAm) , (2-20) 

p=0 
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Figure 3: The diagrams which generate the Majorana couplings. In Pati-Salam language, 
SU(4) x SU{2)l x SU(2)r, the relevant components of the 50(10) representations are: (a) 
il> -> (4, 2, 1), A m -> (TO, 3, 1) as well as (b) ^ -> (4, 1, 2), -> (10, 1, 3). 

with the corresponding diagrams given in Fig. [3j 

While gets a GUT scale VEV in the triplet component of SU(2) R , A^ acquires 
a small induced VEV j2S] 

(Aim) ~ ^ , (2-21) 

in the triplet component of SU(2)l, for details see Appendix|A] The first term in Fig.[3](a) 
then effectively generates a type II see-saw contribution to the physical neutrino mass 
when the Xtb an d the SU(2) L triplet component of A^ get their VEVs. The resulting 
light neutrino mass matrix reads 



^ypeii ~ E<^>'^, (2-22) 

p=0 



with the flavour structure encoded in the three matrices {Xtb)i P = 0; 1> 2. 

On the other hand, the second term of Eq. (12.201) which is depicted in Fig. [3] (6) gives 
rise to the mass of the right-handed neutrinos if A^g acquires a VEV in the direction 

of the right-handed triplet. In addition, the flavon Xtb which is required for PSL2(1) x 
£7(1) invariance also needs to get a VEV. The existence of heavy right-handed neutrinos 
unavoidably leads to a type I see-saw contribution to the light neutrino masses. With the 
neutrino Dirac mass matrix of Eq. (I2.16P , we find 



2 <*&><A° >\~ mr T 



\P=0 / 

The resulting effective light neutrino mass matrix then takes the form 

<ff = Kypel + Kypell ■ (2-24) 

In our PSL 2 (7) model of flavour we want the type II contribution to dominate over type I. 
It is clear from Eqs. ( I2.22|2.23p how this can be achieved. With (x$ ] B ) ~ (A^) ~ M 

both contributions would be more or less equally important. Increasing the scale of (xtb) 
slightly, the type II see-saw contribution increases while, at the same time, the type I 
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see-saw gets suppressed. From now on we will assume dominance of the type II over the 
type I see-saw. The type II see-saw mechanism has the additional advantage of avoiding 
the use of operators with Clebsch zeros in the neutrino direction, a requirement which is 
known to be essential in the type I see-saw mechanism as applied to these models in order 
to provide the necessary suppression required for TB mixing [TT] . 



3 Type II neutrino phenomenology 



In order to extract the neutrino mass spectrum that arises from the first term of Eq. (12.201) , 
we need to insert the VEVs of the flavon sextets Xtb> V = 0, 1, 2 given in Eqs. (l2.7H2.9p 
into Eq. (12.41) . Each of the three flavons comes with its own coupling coefficient so 
that we have to diagonalise 



m 



cff 
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M 2 



p=0 



(3.25) 



As each term is of tri-bimaximal form, we can diagonalise them individually using the 
tri-bimaximal mixing matrix 
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TB 
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(3.26) 
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U T b Ks Utb 



(i-l)Sr < 



2v/3 M 2 
x{ C [°]|(x?L)|-Diag(l,l,-l) 



x 



+ C W|( X ^)|-Diag(l,e^,-e-^) 

+ C [ 2 ]|(x?L)|-Diag(e-^,e^,-l)} , (3.27) 

where the phase factor e ltp is fixed by the PSX 2 (7) group specific parameters 67 and 67 as 

(-3 + iV7) 



br 
b 7 



4 



(3.28) 



numerically corresponding to a phase of (p w 138.6°. Then, after dropping a global phase, 
the three (complex) light neutrino masses rrii are each calculated as the sum of three mass 
parameters 



m 



Ip] 



_ . . C M |/Y bl ^1 



(3.29) 



1 We implicitly absorb the potentially non-zero overall phases of the matrices (xtb) m *° a redefinition 
of the couplings . 



9 



normal neutrino mass ordering 



0.4 
0.35 
0.3 
^0.25 

QJ 

" 0.2 
| 0.15 
. 1 
0.05 




0.05 . 1 0.15 . 2 0.25 . 3 0.35 0.4 

m ml „ (eV) 



inverted neutrino mass ordering 



0.4 
0.35 
0.3 

£■ . 25 

CD 

" 0.2 
| 0.15 
0.1 
. 05 




0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 

mmin (eV) 



Fi gure 4: m ee plotted against m m \ a for the case of a normal mass ordering (left panel, in red) 
and an inverted mass ordering (right panel, in green) assuming for simplicity real parameters 



using different combinations of relative phase factors 
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(3.30) 
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In the special case where the parameters m' p ' are all real, it is possible to find a very 
simple relation with the solar and the atmospheric neutrino mass squared differences, 
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(3.31) 
(3.32) 



As Am^j and |Amjj| tm | have been measured within certain error bars, it is straightforward 
to express m^ 1 ' and as functions of m'°l Note that the sign ambiguity in Am^ tm 
leads to a case distinction between normal and inverted neutrino mass ordering. With 
the (complex) masses rrii taken from Eq. (13.301) we can calculate the effective neutrino 
mass relevant for neutrinoless double beta decay, 



2 1 

- mi + - m 2 



In the case of real m^ they will depend on m' 1 
the smallest (real-valued) neutrino mass 



(3.33) 

only. Furthermore we can easily determine 



m r 



|TOi| 

\m 3 \ 



for normal , 
for inverted 



(3.34) 



neutrino mass ordering which, in the case of real m^, is again only a function of m^. 
We can thus plot m ee against m min , as shown in Fig. |U 

Switching on the complex phases of the parameters m' p l, the dependence of m ee on 
^min becomes much more fuzzy. However, in order to illustrate the phenomenological 
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Figure 5: Scatter plots for m ee (upper panels) and Yli \rrii\ (lower panels) against m m ; n for the 
realistic case of complex parameters for the case of a normal mass ordering (left panels, in 
red) and an inverted mass ordering (right panels, in green). 

consequences arising from the general structure of Eq. (13.301) we have performed a scan 
over the parameters rrM = |mW|e 41/ ' lpl which we have taken to be within the interval 
OeV < |m' p l| < 0.5 eV with arbitrary phases (pW. We have used equal distribution in both 
the absolute values of the masses as well as the phases. With such randomly generated 
input parameters it is straightforward to calculate the resulting complex masses rrti which 
in turn can be converted into atmospheric and solar mass squared differences. Keeping 
only those sets of parameters which lie within the 3a intervals [39] 

< 2.75 x 10~ 3 eV 2 , (3.35) 
7.05xlO~ 5 eV 2 < Am 2 ol = \m 2 \ 2 - \ mi \ 2 < 8.34 x 10~ 5 eV 2 , (3.36) 

we have calculated m ee and m min in order to generate the scatter plot version of Fig. HI 
The result is shown in Fig. We see that allowing for complex phase factors fills in the 
gaps between the branches depicted in Fig. HJ In addition we also show the dependence 
on m m i n of the sum of all light neutrino masses Yli \ m i\ which is relevant for cosmology, 
where the current cosmo logical limit is about Yli \ m i\ < 1-0 e V [40J. 

The histograms in Fig. [6] show the distribution of the numbers of points as a function 
of m min (upper panels) and as a function of m ee (lower panels) corresponding to the scatter 
plots in Fig. \5\ for the case of a normal mass ordering (left panels, in red) and an inverted 
mass ordering (right panels, in green). All the distributions exhibit a very broad peak at 
about ~ 0.05 eV with significant tails out to about ~ 0.4 eV in both m min and m ee . 



2.07 x 10 _3 eV 2 < I Am 



atm I 



"7.3 



mi 
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Figure 6: These histograms show the distribution of the numbers of points as a function of 
mmin (upper panels) and as a function of m ee (lower panels) corresponding to the scatter plots 
in Fig. [5l for the case of a normal mass ordering (left panels, in red) and an inverted mass 
ordering (right panels, in green). These histograms all exhibit a very broad peak at about 
0.05 eV with significant tails out to about 0.4 eV. 



4 Vacuum Alignment in the PSL 2 (7) X 5O(10) model 

We now return to the question of vacuum alignment in the P£X 2 (7) x 5O(10) model, 
using the D-term method and starting with the alignment of the flavon sextets. It turns 
out that, for PSL 2 (7) triplets and anti-triplets, the D-term vacuum alignment method 
reviewed recently in [5] does not work since the invariants involving triplets are the same 
as in £77(3) and the new invariants required for triplet flavon vacuum alignment are not 
present. However, as we shall discuss below, there are new invariants involving sextets 
so the TJ-term method of vacuum alignment is well suited to aligning the flavon sextets. 
Once the sextet flavons are properly aligned, the anti-triplet flavons may then be aligned 
by coupling them to the sextet flavons, as we shall also show. This implies that PSL 2 (7) 
sextets play the crucial role in vacuum alignment quite independently of the crucial role 
that they play in generating third family Yukawa couplings and Majorana masses. 



4.1 Invariants with sextets 

We determine the quartic invariants for the sextet flavon \- Although the sextet irrep 
of PSL>2(7) is real, the physical field need not be real. In general we therefore consider 
invariants of type xlXjXkXi, where we shall assume that such terms involve only sextet 
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flavons of a particular type, i.e. only purely Xtb (^ or a particular choice of p) or the 
flavon Xtop but not operators involving a mixture of different flavon types. This may be 
enforced by assuming a particular messenger sector and introducing a U(l)' symmetry as 
discussed in Appendix [Bj Since we are interested in operators involving only one flavon 
type, we have to consider only the symmetric combinations 

(6©6) s = 1 © 6 © 6 © 8 . 

For x given in the real sextet basis of [2], the representations on the right-hand side can 
be written as 



T 



^ . XiXi 



(4.1) 



/ 



6 : 6 



-2V21xiX5 - QXiXe 
2V / 2Tx2X5 - 6x2X6 
2a/14x 3 X4 + 8x3X6 

- VUxl + 8x4X6 
21x1 + V21 X l ~ 6X5X6 
V-3X? - 3x^ + 4x1 + 4x1 - 3x1 + x\) 



( 



(4.2) 



6 : 0' 



(4.3) 



8 



n 



21x2x3 + 2^x1x4 - 2^x1X6 
21x1x3 + 2^x2x4 - 2^x2X6 

21X1X2 + 2^X3X4 - 2^X3X6 

V7 X i + V7x 2 2 + V7 X i - V7 X l - 2V7 X i - 2V2 X 4X6 
-4^X4X5 - 2^X5X6 
V -V2 X i - V2 X 2 2 - V2 X l - s/2xl - V2xl + ^xl I 

( v / 3X2X3 + X1X4 - 2^X1X5 + 2v / 14XiX6 \ 

v / 2lx2X3 - 3^X1X4 
-VQxl + VQxt - 2X4X5 + 2^X5X6 
-V2~xl - V2xi + 2V2 X l - 2V2~xl + 2^X5 - 2^X4X6 
-v^XiXs + 3^X2X4 
v^XiXs + X2X4 + 2^/6x2X5 + 2a/14X2X6 
-2V2l X 3X5 

\ 2^X1X2-4^X3X4 + 2^X3X6 j 

It should be mentioned that the two sextets and 6' are not defined uniquely since any 
linear combination of G and 0' transforms as a 6 as well. This ambiguity doesn't affect 
the set of quartic invariants even though the explicit form of a particular invariant may be 



(4.4) 
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different. For the octet f2 we have chosen a basis where the generators S^- 8 \ T^ 8 \U^- 8 \ V® 
are real with and diagonal. They can be found in Appendix O The derived 
quartic invariant is necessarily independent of the choice for the octet basis. 

The Kronecker product for the symmetric combinations (6 ® 6) s shows that there are 
six independent quartic invariants of type XiXjXkXi- They can be easily obtained from 
T, 6, 6', Q and their complex conjugates, denoted by Y, 9, 0', Q. Since we have chosen 
real bases, the quartic invariants are formed trivially: 

6 6 

x x = tt, x 2 = ^Te^, x 3 = ^e:e;, (4.5) 

i=l i=l 
6 



* 8=1 i=l 



(4.6) 



x 6 = J2 n - n - • ( 4 - 7 ) 



a=l 



Here we have chosen a convention in which all invariants are real. The way we constructed 
these six invariants obscures a trivial one, namely XiXi) ■ ft is related to the above 
invariants by 



2 



X ° = \ J2xUij = — ^ (49X 1 + 5X 2 + 5X 3 -X 4 + 7X 6 ) . (4.8) 

It is therefore possible to replace Xq in favour of Xo in our set of independent quartic 
invariant s|§ 

4.2 A potential for obtaining sextet flavon alignments 

Let us study a sextet potential of the form [3] 

5 

V = —m 2 ■ a/Xo + A ■ X + A • K a X Q 

a=l 

= -m 2 - y% + (A-Xo)-/, (4.9) 

where we have defined 

f=(l+ E^^V (4-10) 



i X 



In order for this potential to have a minimum, A • / must be positive. Moreover, the 
factor / is independent of the normalisation of (%) so that / is extremised by an ap- 
propriate choice of the alignment. On the other hand, Xq is independent of a particular 



2 In SU(3) the product (6 <E> 6) s ® (6 ® 6) s yields only two independent invariants: To and the sum 
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alignment. Therefore, denoting the extremum of / by fo, the overall scale of the minimum 
is determined from 

](xl)(Xi) = ^/W) 

i=i 



W/vw/ v v-w 2A .^ 



In the following we will assume A > 0. So we need to find the minimum /o > of /. 
However, since we already know the desired alignment vectors, our procedure will be to 
insert these into / and find suitable coefficients n a such that / becomes a (local) minimum. 

4.2.1 Xtb alignment 

Consider the alignment of Eq. ( 12 .7p which breaks PSX 2 (7) down to £4 

(X) oc (0,0,0,0,0,1) T . (4.11) 
We find vanishing first derivatives for each individual invariant, 
d{X a /X Q ) d(l a /l ) 

^75 V" = 7 = °> Va = l,...,5, Vi = l,...,6. 

d{ReXi) d(lmxi) 

In order to get a minimum, the 12 x 12 matrix of second derivatives (the Hessian H a ) 
needs to be positive-definite except for those two real directions which give zero since they 
correspond to the invariance of / under the scaling x% ~~ >* s ' Xit with s 6 C. With the 
alignment (x), this is the Xe direction. For the remaining 10 real dimensions, we obtain 
the following matrices h a 

h x = -8-Diag(0, 0, 0, 0, 0, 1, 1, 1, 1, 1) , (4.12) 

h 2 = 4- Diag (14, 14, 35, 35, 14, 20, 20, 27, 27, 20) , (4.13) 

h 3 = -16 -Diag (14, 14, 14, 14, 14, 9, 9, 9, 9, 9) , (4.14) 

h 4 = -4 -Diag (14, 14, 7, 7, 14, -18, -18, 45, 45, -18) , (4.15) 

K = , k = 28-Diag(2,2,-3,-3,2). (4.16) 

The condition for a minimum of the potential is that all the eigenvalues of the 10 x 10 
matrix 

5 

h = ^ k q h a , 

a=l 

should be positive-definite for appropriate K a . The explicit forms of h a show that many 
potentials can be constructed that are minimised by the alignment of Eq. (14. lip . For 
instance, one could choose 

— 2 < K\ < 10 , K2 — 1 , K3 = K4 = k 5 = . 

The upper bound on k% arises because the sum K\h\ + K 2 /i2 must be positive-definite. 
The lower bound, on the other hand, is related to the requirement of / > 0. Another 
possibility would be 

1 

-— < «3 < , Ki = K 2 = K 4 = K 5 = , 

50 
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where the lower bound arises due to / > 0. We emphasise here that both examples do not 
rely on any tuning and the vanishing of the k coefficients, as assumed in these examples, 
is not a requirement to get a positive-definite h. It is not our intention to give the most 
general potential that is minimised by the alignment (x), but merely to show that a 
potential leading to such an alignment is quite plausible. These examples are sufficient 
to show this. 

Due to the symmetry of such a potential under PSL 2 (7), there exists a discrete de- 
generacy of minima. Given a potential V that is minimised by (x), also the alignment 
Q ■ (x), with Q denoting any element of PSL 2 {7) in the sextet representation, yields a 
minimum. In addition to (%) = (%^), we find six new alignment vectors: 



(X ll] ) oc J -(0,0,0,-714, -721, -1) T , (4.17) 
o 

(X [2] ) oc J -(0,0,0,-714, V2l,-L) T , (4.18) 
o 

(X [3] ) oc -L.(V2l, 721, V2l, V7,0,-V2) T , (4.19) 

(X [4] ) oc ^-(721,-721,-721, V7,0,-V2) T , (4.20) 

(X [5] ) oc -L.(-y/21, V21, -721, 77,0, -72) r , (4.21) 

(X [6] ) oc -i=- (-721, -721,721, 77, 0,-72) T . (4.22) 



The alignment vectors (7 P ') with p = 0, 1, 2 have the form of Eq. (12. 6 p and correspond to 
the vectors (xtb) °f ■^ < l s - (I2.71l2.9p . Therefore, three of the seven discrete minima lead to 
a tri-bimaximal structure. 

It is worth mentioning that the above vacuum alignments (x' p ') can alternatively also 
be obtained form an F-term alignment mechanism. Consider a superpotential of the form 

W = M Xo X + 9X {XX) 6 + g'XoixxYe > ( 4 -23) 

where Xo ^ s a driving sextet field. The parentheses denote the contraction to the two 
distinct sextet as shown in Eqs. ( I4.2|4.3p . It is easy to show that the resulting F-term 
equations of Xo are solved by the sextet alignment (x' ') = (0, 0, 0, 0, 0, v) T if 

M 

9 + 572> ' 

Due to the symmetry of the superpotential under PSL 2 (7), the other six alignments 
(x^) lead to vanishing F-terms as well. Thus the three sextets leading to a tri-bimaximal 
structure can be aligned alternatively through an F-term mechanism. In the case of 
charged sextets the mass parameter Ai must be generated dynamically by the VEV of 
some PSL 2 (7) singlet. 
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4 - 2 - 2 Xtop alignment 

We now turn to the study of the alignment of Eq. ( 12.101) . In general we seek to solve 
ten non-trivial conditions with a set of five parameters K{. Therefore it is by no means 
guaranteed that an extended parameter space of solutions to the minimisation conditions 
of the potential should exist. Remarkably, by plugging Eq. (I2.10p into the potential 
in Eq. (14. 9p . we find that, for a certain choice of parameters discussed below, the first 
derivatives of / can be made to vanish. However, unlike the previous alignments, first 
derivatives of / do not vanish for all the invariants taken separately. 

To be precise, a straightforward calculation shows that, in order to achieve the desired 
top alignment (Xtop)j the vanishing of the first derivatives requires only two relations 
amongst the different sextet combinations which result from the symmetric combinations 
of two sextets, to wit n 2 = k 3 = k 4 + n 5 /\^f. Clearly the origin of these relations remains 
to be understood, and, for example, could result from some underlying higher symmetry, 
although this goes beyond the scope of the present PSL 2 (7) discussion. However we 
emphasise that the fact that a successful potential can be found with particular values of 
parameters which can lead to the desired top alignment is highly non-trivial and this is 
not the case in general for other alignmentso 

Assuming the above two relations, the potential that has a chance of being minimised 
by (Xtop) includes an / of the form 



X + k X x + k' (X 2 + X 3 + X 4 ) + k" fx 4 - V7Z, 
^-v"" «■ „ ' A 



= X =X' ^~ 

=X" 



(4.24) 



Since / has to be positive, we get the condition 
< / 



70 140 

= 1 + • k + ^ • k' + • k" . (4.25) 

(xtop) a 



On the other hand, the Hessian needs to be positive-definite. The calculation for all three 
independent contributions to / shows that the 10-dimensional Hessian h has two positive 
and eight zero eigenvalues, while h! has six positive and four negative ones, and finally h" 
has ten negative eigenvalues. We therefore conclude that the choice 

k = k' = , -A < k" < , (4.26) 

leads to a potential V which is minimised by the alignment of Eq. ( I2.10p . In other words, 
X" must be suppressed compared to X . Other solutions are possible as well, but as the 
three Hessians cannot be diagonalised simultaneously, it is not easy to give a general 
expression. However we emphasise that K and k' need not be zero, for instance the choice 

K = K> = -«"=_, 



3 To illustrate this, consider the alignment (x a it) ^ ( — b 0' 0' V3/b7, \/2, 0) T which leads to a Yukawa 
matrix with identical (2, 2) and (3, 3) elements and zero entries everywhere else. Requiring vanishing hrst 
derivatives fixes the parameters Ki uniquely. Up to an overall scale we find n\ — 0, k-i — K3 = 2^4 = 
2kc,/\/7. However, in this case, the matrix of second derivatives turns out to have positive and negative 
eigenvalues showing that no quartic potential exists which is minimised by (x a it) - 
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leads to an acceptable potential as well. 

Having illustrated that the alignment (x top ) can be obtained from reasonable poten- 
tials, we again find that a transformation of this alignment vector under PSL 2 {7) yields 
new minima. In this case, we get 168 different minima, including the one of Eq. (12.101) . 
All of them are physically identical since we can use the freedom to redefine our basis by 
applying a suitable PSL 2 (7) symmetry transformation. 

4.3 Flavon anti-triplet vacuum alignment 

As outlined in [2] we need to introduce anti-triplet flavon fields to generate the first and 
second family Yukawa couplings. The question arises how to align such anti-triplets of 
PSL 2 {7). The immediate idea to study quartic terms of the form <p\<pi <p\(f)\ which proves 
successful for flavour groups like A 2 7 [H], Z 7 x Z 3 [20], A 4 [7] and others [3J, does not lead 
to an alignment of the anti-triplets because PSL 2 (7) is too big a symmetry so that the 
trivial SU(3) invariant 0j'0 l 0^0 : ' is the only allowed quartic term of only anti-triplets. 

However, due to the presence of the sextet representation, the possibility of aligning 
the anti-triplets by coupling them to the pre-aligned sextets arises. It turns out that 
the simplest operator Xtb $ * s n °t suitable. After inserting the VEV of Xtb we obtain 
a contribution to the anti-triplet scalar potential which is extremised by the alignment 
(0) oc (1, x, x). Due to the degeneracy of vacua the desired anti-triplet alignments of 
Eq. (12. lip cannot be generated form terms like 0XtL ^ where p = 0, 6@ 

Since the simplest operator discussed above does not work, in the following we discuss 
the next simplest possibility of coupling two anti-triplets to the product of two sextet 
flavons that have the tri-bimaximal alignment, (Xtb) wr th p = 0, 1,2. We then prevent 
the flavon anti-triplets from coupling to the sextet flavon x top by means of a particular 
messenger sector and a U(l) symmetry as discussed in Appendix [B] The relevant scalar 
potential for the anti-triplets then reads 

V = -m 2 .j2<PW + \(i2^A +tEE c ^^*l(xS) i . (4-1) 

i=l \i=l / a=l i,j=l k,l=l 

Here the index a labels the three different invariants of the PSL 2 {7) product 

(3 g> 3) <g> (6 <g> 6) , 

108 l©2-6©7e2-8 

with the index structure of c a '* fcZ being defined by the Clebsch Gordan coefficients of the 
corresponding Kronecker products. The contraction to the singlet yields a term of the 
form 

3 

AV = a -J2& , ( 4 -2) 

i=l 

4 Note, however, that it is possible to use similar couplings to align anti-triplets against the T preserving 
sextet Xt w ith an alignment of the form (\t) °z (v^, 0, 0, 0, 1, z) T , where z remains undetermined. It 
turns out that the resulting anti-triplet alignment (</>) oc (1,0,0) preserves T as well. 
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where ao includes the VEVs of the sextet flavons. This contribution to the potential 
does not constrain the anti-triplet alignment. The situation changes for the other two 
invariants obtained form the contractions to the octet. Inserting the alignments (x*rL) 
and (xtb)> whh p,p' = 0, 1, 2, into Eqs. (1C.5|IC6|) we find that, in general, only the third 
and the forth components of the octet of 3 <8> 3 in Eq. (1C7[) can be projected out in the 
potential. However, with p = p', i.e. identical sextets entering in the potential, both terms 
vanish identically. Only in the case where p ^ p 1 we get non-zero contributions to the 
potential. Choosing for example p — 1 and p' — 2, the symmetric octet of Eq. (1C.5I) is pro- 
portional to (0, 0, 0, 1, 0, 0, 0, 0) while the antisymmetric octet of Eq. (IC.6I) is proportional 
to (0, 0, 1, 0, 0, 0, 0, 0). We thus find two additional independent terms 



AV a = a a 



(4-3) 

I(-20 1 + 2 + 3 ) + ^(0 1 + 2 -20 3 ) + ^(0i_ 2 2 + 3 )l . (4.4) 



Combining these two terms linearly with each other and with AVo we can define two new 
independent terms which determine the anti-triplet alignment^! 

AVi = « 1 (0 t 1 + ^ + ^)(0l + 02 + ^3 ) ? (45) 

AV 2 = a 2 (^-^)(0 2 -0 3 ) . (4.6) 

These two terms of the scalar potential are at the core of the discussion of the anti- 
triplet alignment. Choosing the values of a± and a 2 appropriately gives rise to a potential 
which is minimised by an alignment of the anti-triplets required to generate the first and 
second family Yukawa couplings. As an aside we note that the two independent terms of 
Eqs. ( I4.5|4.6p can be obtained similarly from Eq. (14. ip using tri-bimaximal flavon sextets 
with p = and p' = 1 or 2. 

Let us first consider the two corrections to the potential individually. For positive a% 
the minimum of AVi is zero. This entails a partial alignment of the form 

x 

ai > : ((f)) oc | y 

-x - y / 

where x, y G C remain undetermined. In contrast, for negative ol\ the resulting alignment 
is completely fixed 

- ( v 

oci < : ((f)) oc 1 

Similarly, the potential term AV2 gives rise to the following structure of the anti-triplet 
VEVs 

/x\ /0 

a 2 > : (<p) oc \y , a 2 < : ((f)) oc 1 

W v-i 



5 Explicitly, AVy/ai = AV s /a s + AV /a Q and AV 2 /a 2 = \{AV a /a a - AV s /a s + 2AV /a ). 
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Remarkably, the alignments derived from both terms AVi and AV 2 individually can be 
made compatible by choosing the signs of oi\ and a 2 according to the following combina- 
tions 



Thus it is possible to have an anti-triplet flavon field 023 whose VEV becomes aligned 
along (0, 1, —1), while another flavon field 0i 2 3 ends up with the alignment (1, 1, 1). 

5 Conclusion 

In this paper we have constructed a realistic SUSY GUT of Flavour based on PSL 2 (7) x 
5*0(10), where the quarks and leptons in the 16 of 5*0 (10) are assigned to the complex 
triplet representation of PSL 2 (7), while the flavons are assigned to a combination of 
sextets and anti-triplets of PSL 2 (7). It represents the first model based on the finite 
group PSL 2 {7), which is the smallest simple group that contains both complex triplet 
and (real) sextet representations. This group seems particularly well suited to 50(10) 
since the sextets may be used to provide the large third family Yukawa coupling, as well 
as type II neutrino masses. Furthermore PSL 2 {7) contains 5 4 as a subgroup, and this 
allows the possibility of explaining TB neutrino mixing in a direct way, by preserving 
the generators S,U of PSL 2 (7) in the neutrino sector, which become identified as the 
neutrino flavour symmetry. There are very few models that can account for TB neutrino 
mixing using 5O(10) and this is the first model which can do this directly. 

Using a D-term vacuum alignment mechanism, we have shown how the flavon sextets 
of PSL 2 (7) can be aligned along the 3-3 direction leading to the third family Yukawa 
couplings. Such sextets aligned along the 3-3 direction are ideally suited for giving a 
universal contribution to the 3-3 Yukawa coupling at the lowest possible one-flavon order 
in 5O(10) models, allowing a sizeable universal top-bottom-tau Yukawa coupling. We 
emphasise that the fact that a successful potential can be found with particular values 
of parameters which can lead to the desired 3-3 vacuum alignment of the sextet flavons 
is highly non-trivial and this is not the case in general for other alignments. However, in 
order to realise the flavon sextet potential that yields an alignment along this 3-3 direction, 
it is necessary to assume certain relations amongst the parameters of the potential. These 
relations could in principle emerge from a higher symmetry, beyond PSL 2 (7), although 
this takes us beyond the scope of the present paper, though it should be the subject of 
future investigation. 

Other sextets are aligned along the neutrino flavour symmetry preserving directions 
in an even more natural way, without requiring any relations between the parameters of 
the potential, and such alignments suggest the possibility of TB neutrino mixing via a 
type II see-saw mechanism. We have explored the phenomenological consequences of such 
a type II see-saw mechanism and obtained statistical predictions for neutrinoless double 
beta decay and neutrino masses in cosmology. The distributions of randomly generated 
points exhibit a very broad peak at about ~ 0.05 eV with significant tails out to about 
~ 0.4 eV in both m min and m ee . 




(4.7) 
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Anti-triplet flavons are also introduced and aligned against the pre-aligned TB sextet 
flavons, in order to give the remaining structure of the charged fermion mass matrices. In 
principle these may also be used to account for neutrino masses and TB mixing via a type 
I see-saw mechanism as in [TTJ, but since this is well known we have focused on the new 
type II possibility opened up by having sextets in the model. This also avoids the use of 
operators with zero Clebsch factors in order to provide the necessary suppression in the 
type I Dirac neutrino sector [TTJ. Nevertheless, the anti-triplet flavons are instrumental 
in giving the successful mass matrices in Eq. (12.181) via the assumption of two different 
messenger mass scales in the up and down sectors. In this model the mass matrices are 
achieved in a very natural way since, with the inclusion of the singlet flavon £, the first 
row and column is cubic in the messenger mass, while the 2-3 block is quadratic and the 
3-3 element involves the universal Higgs messenger mass Ma- 
in conclusion the SUSY GUT of Flavour based on PSL 2 (7) x 50(10) leads to a very 
elegant model, combining the mathematical beauty of PSL2(7) with the attractiveness 
of SO (10) unification, and solving many of the problems related to achieving successful 
fermion masses and TB mixing in the SO (10) framework. The SUSY flavour problem 
is also solved here exactly as in the SU(3) model discussion in [41j, since both models 
use identical anti-triplet flavon alignments. Finally we emphasise that the type II see-saw 
mechanism in this model is consistent with neutrinoless double beta decay right up to the 
limit of current experiments. 
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Appendix 



In 50(10) the induced VEV of the Ay^ field in Fig. |3](a) is obtained by replacing the 
Aj26 leg of the diagram by 



A The induced VEV 
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The relevant components of the 5*0(10) representations are 



H w (15,2,2) , 5 126 -+ (10,1,3) , (A.l) 

with if 10 acquiring a VEV at the electroweak scale while S 12 6 gets a GUT scale VEV. 
Assuming the messengers in the complete diagram to have a mass of order M, the induced 
VEV can be calculated to be 

( S 12e) (10,1,3) (Hio) (1,2,2) (H 10 ) (i )2 ,2) , . 

~M~ 2 M ' (A ' 2) 



where, for simplicity, we have taken (Hi2e)(ioi3) ~ It is clear from Eq. (1A.2|) that 



the resulting contribution to the neutrino masses corresponds to the type II see-saw 
mechanism because the two left-handed doublets of the H w pair combine to an SU(2)l 
triplet, whereas the two right-handed doublets are contracted with the right-handed triplet 
of S126 to a singlet under SU(2)r. The relevant component of the messenger is thus 
an 577(2) j, triplet. 



B Controlling the flavon potential 

In Section S] we discussed the terms of the flavon potential that are required to obtain the 
alignments that generate the Majorana and Yukawa couplings. There we encountered the 
following four types of couplings 



t t [pit [p] [pit [p] 

Xtop Xtop Xtop Xtop > Xtb Xt bXtbXtb > 



^123^123XtbXtB ' 



(B.l) 



which involve four sextet and two anti-triplet flavon fields. As we have already discussed 
the contractions that yield P5L 2 (7) invariants, we suppress all indices in Eq. (IB.ip . At 
this point it is important to observe that the suggested vacuum alignment is based on the 
absence of terms like 



t [p] t [p] 

Xtop Xtop Xtb Xtb ' 



XtB^XtbXtB^XtB > 023 < / , 23XtopXt' 



op A top ) 



0123 < / ) 123XtopXtop > 

(B.2) 



that - at the effective level - cannot be forbidden by symmetries alone. It is therefore 
necessary to resort to a particular messenger sector. A neutral messenger would automat- 
ically give rise to diagrams such as 



0t 



which cannot distinguish between the structure of the terms in Eq. fIB.ip and Eq. (IB. 21) . 
Therefore the messengers must be charged under additional symmetries. Accordingly, a 
possible way to forbid the terms of Eq. (1B.2j) is given by diagrams of type 



22 



where the charge of the respective messengers determines whether or not x' — X is allowed 
as well as which sextet flavons can couple to the anti-triplet flavons. Most notably we 
need to separate the top sextet from the sextets of tri-bimaximal type. This is achieved 
easily by assigning different U(l) charges q to the flavons Xto P anc l Xtb- Choosing for 
instance 

q( Xtop ) = -1 , q( X ? B ) = "2 , <?(0 23 ) = "2 , <?(0 123 ) = 4 , (B.3) 

one can generate the operators of Eq. fIB.ip using messengers with charges q — 2,4, 4, —2, 
respectively, while the first, third and fourth term of Eq. (IB. 21) would require messengers 
with odd U(l) charge. In the absence of such messengers the top sextet flavon cannot 
mix with the tri-bimaximal one. 

Furthermore, we also need to forbid the second term of Eq. (1B.2|) . i.e. the mixing 



among the three tri-bimaximal flavon sextets. For this purpose we introduce a separate 
symmetry, U(l)', which distinguishes between Xtbi with P — 0,1,2. One possible set of 
U(l)' charges could be 

qXxfs) = 2 , q\x [ T ] B) = 3 , JiA) = 5 • (B.4) 

The second term of Eq. (IB.IJ . corresponding to three distinct quartic operators, would 
arise from messengers with charges q' = —4, —6, —10, respectively. At the same time, the 
analogous mixing term of Eq. flB.2j) would require messengers with either q' = —5, —7, —8 



or q' = 1,2,3. 

In the construction of the complete model it is necessary to introduce three PSL2(7) 
singlet flavons These are associated with the three tri-bimaximal sextet flavons Xtb 
and carry opposite U(l)' charge. Then the Xtb l e S s °f the diagrams in Fig. [3] need to be 
replaced by 



Wb 



With this completion, the only particles that are charged under the U(l)' symmetry are 
Xtbi C' p ' an d the above mentioned messengers with q' = —4, —6, —10. 



C The octet representation 

Just like in SU(3), the octet can be constructed from the product 3 ® 3, see for instance 
Ref. [31 J. The resulting 8x8 matrices for the PSL2(7) generators are real, but none 



23 



is diagonal. Since we are interested in the combination (6 <S> 6) — > 8, it is convenient to 
perform a similarity transformation to a basis in which - analogous to the sextet (see [2]) - 
the generators and W' 8 ' are diagonal. We obtain 
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(C.2) 



(C.3) 



(C.4) 



With the identification 

A = [r [8] W [8l 5 [8l (T [8] ) 2 W [8] ] _1 V [8] [T [8] U [8] S [8] (T [8] ) 2 l{ [8] ] 



B 



one can easily check that the presentation of PSL 2 (7) is satisfied. The symmetric octet Q 
of Eq. f)4.4p is given in the basis of Eqs. flC.lHC.4]) . More generally there are two inde- 
pendent octets derived from the product of two sextets which we take in the basis of [2]. 
They are 



/ VS X2 X3 + X1X4 - 2v%iX5 + 2VU XlX ' 6 

V2IX2X3 -3V7x x X4 
-VGXiX'i + V6X2X2 - 2X4X5 + 2^14X5X6 
V2XxXi " V2X2X2 + 2^X3X3 " 2^X4X4 + 2^X5X5 

-V2IX1X3 + 3^X2X4 
v^XiXs + X2X4 + 2^X2X5 + 2^X2X6 
-2V2IX3X5 

\ 2^X1X2 - 4^X3X4 + 2^X3X6 



\ 



2%/7x 4 X6 



+ (x x'), 



(C.5) 
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f-V2Tx 2 x' 3 - VIXiX'a ~ 2v / 2XiX 6 \ 
V3x 2 x' 3 ~ 3 XiX 4 - 2V6X1X5 
-2VW5 - Sv^XsXe 
-6X4X6 

-v^XiXs + 3x 2 X 4 - 2^6x2X5 
21X!X3 " V7 X2 X4 ~ 2^X2X6 
2^X1X2 + 2^X3X5 
6X3X6 



- (x x') 



(C.6) 



Obviously, the antisymmetric one only exists for x 7^ Adopting the triplet basis of [2], 
the octet of the product 3 <g> 3 reads 



/ fa (40 1 + 2 + 3 ) + (fa + 3 ) (0 1 - 20 2 - 20 3 ) \ 

-3i[fa(fa + fa)-(fa + fa)fa] 
-v / 2[0 1 (-20 1 + 2 + 3 ) + fa(fa + fa- 20 3 ) + 03(0 X - 20 2 + fa)} 

a/6[0i(0 2 + P) + MP + P) + MP + P)\ 

V3[fa(P - ^ 3 ) + + 20 2 ) - 03(0' + 20 3 )] 

iV3[^i(^ - 3 ) - fa(P + 20 3 ) + 3 (0 X + 20 2 )] 

VeifaiP - P) + MP - 2 ) - MP - P)} 



\ 
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